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Abstract. This paper concerns the reconstruction of an anisotropic conduc- 
tivity tensor in an elliptic second-order equation from knowledge of the so-called 
power density functionals. This problem finds applications in several coupled- 
physics medical imaging modalities such as ultrasound modulated electrical 
impedance tomography and impedance-acoustic tomography. 

We consider the linearization of the nonlinear hybrid inverse problem. We 
find sufficient conditions for the linearized problem, a system of partial differen- 
tial equations, to be elliptic and for the system to be injective. Such conditions 
are found to hold for a lesser number of measurements than those required in 
recently established explicit reconstruction procedures for the nonlinear prob- 
lem. 

1. Introduction. In the context of hybrid medical imaging methods, a physical 
coupling between a high-contrast modality (e.g. Electrical Impedance Tomography, 
Optical Tomography) and a high-resolution modality (e.g. acoustic waves, Mag- 
netic Resonance Imaging) is used in order to benefit from the advantages of both. 
Without this coupling, the high-contrast modality, usually modeled by an inverse 
problem involving the reconstruction of the constitutive parameter of an elliptic 
PDE from knowledge of boundary functionals, results in a mathematically severely 
ill-posed problem and suffers from poor resolution. The analysis of this coupling 
usually involves a two-step inversion procedure where the high-resolution modal- 
ity provides internal functionals, from which we reconstruct the parameters of the 
elliptic equation, thus leading to improved resolution [2, 5, 14, 17, 22]. 

A problem that has received a lot of attention recently concerns the reconstruc- 
tion of the conductivity tensor 7 in the elliptic equation 

V-(7Vu)-0 {X), u\ax=9, (1) 
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from knowledge of internal power density measurements of the form Vu-^Vv, where 
u and V both solve (1) with possibly different boundary conditions. This problem 
is motivated by a coupling between electrical impedance imaging and ultrasound 
imaging and also finds applications in thermo-acoustic imaging. 

Explicit reconstruction procedures for the above non-linear problem have been 
established in [9, 6, 20, 19, 17], successively in the 2D, 3D, and nD isotropic case, 
and then in the 2D and riD anisotropic case. In these articles, the number of func- 
tionals may be quite large. The analyses in [17] were recently summarized and 
pushed further in [18]. If one decomposes 7 into the product of a scalar function 
T = (det7)" and a scaled anisotropic structure 7 such that det7 = 1, the lat- 
ter reference establishes explicit reconstruction formulas for both quantities with 
Lipschitz stability for t in and involving the loss of one derivative for 7. 

In the isotropic case, several works study the above problem in the presence 
of a lesser number of functionals. The case of one functional is addressed in [4], 
whereas numerical simulations show good results with two functionals in dimension 
n = 2 [1, 12]. Theoretical and numerical analyses of the linearized inverse problem 
are considered in [15, 16]. The stabilizing nature of a class of internal functionals 
containing the power densities is demonstrated in [16] using a micro-local analysis 
of the linearized inverse problem. The above inverse problem is recast as a system 
of nonlinear partial differential equations in [3] and its linearization analyzed by 
means of theories of elliptic systems of equations. It is shown in the latter reference 
that 71-1-1 functionals, where n is spatial dimension, is sufficient to reconstruct a 
scalar coefficient 7 with elliptic regularity, i.e., with no loss of derivatives, from 
power density measurements. This was confirmed by two-dimensional simulations 
in [7]. All known explicit reconstruction procedures require knowledge of a larger 
number of internal functionals. 

In the present work, we study the linearized version of this inverse problem in 
the anisotropic case, i.e. we write an expansion of the form 7' = 70 + £7 with 
70 known and e <C 1, and study the reconstructibility of 7 from linearized power 
densities (LPD). We first proceed by supporting the perturbation 7 away from the 
boundary dX and analyze microlocally the symbol of the linearized functionals, and 
show that, as in [16], a large enough number of functionals allows us to construct a 
left-parametrix and set up a Fredholm inversion. The main difference between the 
isotropic and anisotropic settings is that the anisotropic part of the conductivity is 
reconstructed with a loss of one derivative. Such a loss of a derivative is optimal 
since our estimates are elliptic in nature. It is reminiscent of results obtained for a 
similar problem in [8]. 

Secondly, we show how the explicit inversion approach presented in [17, 18] carries 
through linearization, thus allowing for reconstruction of fully anisotropic tensors 
supported up to the boundary oi X. In this case, we derive reconstruction formulas 
that require a smaller number of power densities than in the non-linear case, giving 
possible room for improvement in the non-linear inversion algorithms. 

For additional information on hybrid inverse problems in other areas of (mostly 
medical) imaging, we refer the reader to, e.g., [2, 5, 21, 22]. 



2. Statement of the main results. Consider the conductivity equation (1), 
where X C K" is open, bounded and connected with n > 2, and where 7"^ is a 
uniformly elliptic conductivity tensor over X. 
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We set boundary conditions {gi, . . . ,gm) and call uf the unique solution to (1) 
with wf lojf = 1 < i < TO and conductivity 7^. We consider the measurement 
functionals 

H!^:j'^H!j{j') = Vul-j'^u'^{x), l<*,j<m, x e X. (2) 

Considering an expansion of the form 7" = 70 + £7, where the background 
conductivity 70 is known, uniformly elliptic and e so small that the total 7'^ remains 
uniformly elliptic, we first look for the Frechet derivative of (2) with respect to 7 
at 7o. Expanding the solutions uf accordingly as 

^ Ui + evi + 0{e'^), 1 < i < m, 

the PDE (1) at orders ©(e") and 0{e^) gives rise to two relations 

-V-(7oVu,) = (X), u^\ax=g^, (3) 

-V • (7oVw.) = V • (7VUO (X), v.lox^O- (4) 
The measurements then look like 

H^j = Vu, ■ 7oVuj + e (Vu^ • jVuj + Vu, ■ 7oVwj + Vuj ■ jqVv,) + ©(e^). (5) 
Therefore, the component dHij of the Frechet derivative of H at 70 is 

dHij{'j) = Vui ■ jVuj + ■ 7oVwj + Vuj • 70 Vw;, x G X, (6) 

where the Wi's are linear functions in 7 according to (4). 

In both subsequent approaches, reconstruction formulas are established under 
the following two assumptions about the behavior of solutions related to the con- 
ductivity of reference 70. The first hypothesis deals with having a basis of gradients 
of solutions of (3) over a certain subset il C X. 

Hypothesis 2.1. For an open set C_ X, there exist (171, . . . , ,g„) £ H^{dX)" 
such that the corresponding solutions (ui, . . . of (3) with boundary condition 
Uildx ^ gi satisfy 

inf dct(Vui, . . . , Vii„) > cq > 0. 

Once Hypothesis 2.1 is satisfied, any additional solution Un+i of (3) gives rise to 
a, n X n matrix 



Z =[Zi\...\Zn], where := V — — — ^ . (7) 

det(Vwi,...,Vu„) 

As seen in [17, 18], such matrices can be computed from the power densities {Vui ■ 
7oV7ij}"^^2^ and help impose orthogonality conditions on the anisotropic part of 
70- Once enough such conditions arc obtained by considering enough additional 
solutions, then the anisotropy is reconstructed explicitly via a generalization of the 
usual cross-product defined in three dimensions. In the linearized setting, wc find 
that one additional solution such that Z has full rank is enough to reconstruct the 
linear perturbation 7. We thus formulate our second crucial assumption here: 

Hypothesis 2.2. Assume that Hypothesis 2.1 holds over some fixed C X. There 
exists gn+i G H^{dX) such that the solution of (3) with boundary condition 
Un+i\dx = 9n+i has a full-rank matrix Z (as defined in (7)) over 57. 
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Remark 1 (Case 70 constant). In the case where 70 is constant, then it is straight- 
forward to sec that gi = Xi\gx (1 < « < 't-) f^ulfill Hypothesis 2.1 over X. Moreover, 
if Q = {Qij}i j=i denotes an invertible constant matrix such that Q : 70 = 0, then 
the boundary condition 9,1+1 := ■^qijXiXj\gx fulfills Hypothesis 2.2, since we have 
Q = Z. 

Throughout the paper, we use for (real-valued) square matrices A and B the 
contraction notation A : B = ti AB^ = J^i j ^ijBij, with B^ the transpose matrix 
of A 

Remark 2. In the treatment of the non-linear case [6, 20, 17, 18], it has been 
pointed out that Hypothesis 2.1 may not be systematically satisfied globally in 
dimension n > 3. A more general hypothesis to consider would come from picking 
a larger family (of cardinality > n) of solutions whose gradients have maximal rank 
throughout X. While this additional technical point would not alter qualitatively 
the present reconstruction algorithms, it would add complexity in notation which 
the authors decided to avoid; see also [8]. 

2.1. Past work and heuristics for the linearization. In the reconstruction 
approach developpcd in [19, 17, 18] for the non-linear problem, it was shown that not 
every part of the conductivity was reconstructed with the same stability. Namely, 
consider the decomposition of the tensor 7' into the product of a scalar function 
T = (det7')" and a scaled anisotropic structure 7' with det7' = 1. The following 
results were then established. Starting from n solutions whose gradients form a 
basis of R" over a subset fi C X, it was shown that under knowledge of a W^'°°{X) 
anisotropic structure 7', the scalar function logdet7' was uniquely and Lipschitz- 
stably reconstructible in W^'°°{il) from W^'°° power densities. Additionally, if one 
added a finite number of solutions Un+i, . . . , Un+i such that the family of matrices 

. . . , defined as in (7) imposed enough orthogonality constraints on 7', then 
the latter was explicitely reconstructible over from the mutual power densities of 
(wi, . . . , Un+i)- The latter reconstruction was stable in L°° for power densities in 
W^'°° norm, thus it involved the loss of one derivative. 

Passing to the linearized setting now (recall 7*^ = 70 + £"f), and anticipating 
that one scalar quantity may be more stably reconstructible than the others, this 
quantity should be the linearized version of logdet7'^. Standard calculations yield 

logdet(7o -t- £7) = logdet7o + logdet(I„ + £7,7^) = logdet7o + etr (7(7 V) + 0{e'^), 

and thus the quantity that should be stably reconstructible is tr (7o^^7). The lin- 
earization of the product decomposition (t, 7') above is now a spherical-deviatoric 
one of the form 

7 = ^tr (7c7'7)7o + 7d, Id ■= loho'^f"^ (») 
where '^'^^ is the linear projection onto the hyperplane of traceless matrices A'^'^^ := 

2.2. Microlocal inversion. The above inverse problem in (4)-(6) may be seen as 
a system of partial differential equations for (7, {wj}). This is the point of view 
considered in [3]. However, {vj} may be calculated from (4) and the expression 
plugged back into (6). This allows us to recast dH as a linear operator for 7, 
which is smaller than the original linear system for (7, {vj}), but which is no longer 
differential and rather pseudo-differential. The objective in this section is to show. 
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following earlier work in the isotropic case in [16], that such an operator is elliptic 
under appropriate conditions. 

Wc first fix il' CC X and assume that supp7 C ft', so that integrals of the 
form Jg„ e''^'^p(x,^) : 7(0'^^ sltc well-defined, with p{x,£^) a matrix-valued symbol 
whose entries are polynomials in ^ (see [11, p. 267]) and where the hat denotes the 
Fourier Transform 7(^) — Jg^„ e~''''^'^j{x) dx. We also assume that 70 G C°°{il') 
and can be extended smoothly by 70 ~ outside fl' . As pointed out in [16], in 
order to treat this problem microlocally, one must introduce cutoff versions of the 
dHij operators, which in turn extend to pseudo-differential operators (^'DO) on K". 
Namely, if ft" is a domain satisfying f2' CC fJ" CC X and xi is a smooth function 
supported in X which is identically equal to 1 on a neighborhood of Q" , the operator 
7 i-> XidHij{xil) can be made a ^fDO upon considering Lq = —V • (70V) as a 
second-order operator on M" and using standard pseudo-differential parametrices to 
invert it [13]. We will therefore not distinguish the operators dHij from their pseudo- 
differential counterparts. The task of this section is then to determine conditions 
under which a given collection of such functions becomes an elliptic operator of 7 
over fi'. 

Using relations (4) and (6), we aim at writing the operator dHij hi the following 
form 

dH,,{x)^{2i,)-- ( [ e'-^-'-^-yHu,{x,0--l{y)didy, (9) 

with symbol My (x,^) (pseudo-differential terminology is recalled in Sec. 3.1). We 
first compute the main terms in the symbol expansion of dHij (call this expansion 
Mij — Mij\Q + A/y l-i + 0(1^1"^) with Mij\p homogeneous of degree p in ^). From 
these expressions, we then directly deduce microlocal properties on the correspond- 
ing operators. 

The first lemma shows that the principal symbols M^ lo can never fully invert for 
7, no matter how many solutions Ui we pick. When Hypothesis 2.1 is satisfied, then 
the characteristic directions of the principal symbols {M^ (x, ^)}i<i j<„ reduce to 
an— 1-dimensional subspace of S'„(K). Here and below, we recall that the colon 
":" denotes the inner product A : B = ti {AB^) for (A, B) G S'„(M) and denotes 
the symmetric outer product U QV = \{U ®V + V ®U) lor U,V . 

Lemma 2.3. (i) For any i,j and x Cz X , the symbol Mij\o satisfies 

Mtjlo ■■ (70? V) = 0, for all ij e S"-^ satisfying ■ ? 0. (10) 

(ii) Suppose that Hypothesis 2.1 holds over some f2 C X. Then for any x ^ Q, if 
P £ 5,1 (M) is such that 

My|o:P = 0, l<i<j<n, (11) 

then P is of the form P = 70? 77 for some vector rj satisfying 77 • ? = 0. 

Since an arbitrary number of zero-th order symbols can never be elliptic with 
respect to 7, we then consider the next term in the symbol expansion of dHij. We 
must also add one solution Un+i to the initial collection, exhibiting appropriate 
behavior, i.e. satisfying Hypothesis 2.2. The collection of functionals we consider 
below is thus of the form 

dH := {dH^j \ l<i<n, i < j <n + l}, (12) 

and emanates from n -|- 1 solutions (wi, . . . ,w„+i) of (3) satisfying Hypotheses 2.1 
and 2.2. 
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In order to formulate the result, we assume to construct a family of unit vector 
fields 

homogeneous of degree zero in ^, smooth in x and everywhere orthonormal. We 
then define the family of scalar elliptic zeroth-order VPDO 

T : J ^ = {Tpqj}o<p<q<n-i, whcrc 
Tpq-fix) (27r)-" / e'-^'^A^'ip i^A^' : 7(C) d^, I < p < q < n, ^^^^ 

which can be thought of as a microlocal change of basis after which the operator 
dH{j) becomes both diagonal and elliptic. Indeed, we verify (see section 3.4) that 
for any fc > 1 and 7 sufficiently regular, we have 

||7||h'=(o') < C||T7||ff.(^,) + C^hWmnn < C^hWHHn')- (14) 

The above estimates come from standard result on pseudo-differential operators 
[13]. The presence of the constant C2 indicates that T can be inverted microlocally, 
but may not injective. 

Composing the measurements dHij with appropriate scalar Vl/DO of order and 
1, we are then able to recover each component of the operator (13). The well- 
chosen "parametrices" are made possible by the fact that the collection of symbols 
Mij\o + Mijl^i becomes elliptic over il' when Hypotheses 2.1 and 2.2 are satisfied. 
Rather than using the full collection of measurements dH (12), we will consider the 
smaller collection {(ii?ij}i<ij<n augmented with the n measurement operators 

n 

Liil) = + ^ dHi^ri+i{'y), 1 < i < n, (15) 

where (/ii, . . . , fin, fj){x), known from the measurements {Hij}^J]^^, are the coeffi- 
cients in the relation of linear dependence 

/iiVui + • ■ • + /-f„Vu„ + /iVu„+i = 0. 

We also define the operator G with principal symbol — i.||^o?||- Our conclu- 
sions may be formulated as follows: 

Proposition 1. Let the measurements dH defined in (12) satisfy Hypotheses 2.1 
and 2.2. 

(i) For (a,/3) = (0,0) and l<a<P<n — \, there exist {Qai3ij}i<i<j<n G 
such that 

Qc.prj°dH,j=Tap mod*-i. (16) 

l<i,^'<n 

(a) For any 1 < a < n — 1, there exist {Bai\i<i<n G '1'° such that the following 
relation holds 

lI oBa;oL,- RaoR^Toa mod*-\ (17) 

where the remainder Ra o R can be expressed as a zeroth-order linear combi- 
nation of the components Too cind {Tpq}Kp<q<n-i reconstructed in (i). 
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The presence of the Lg term in part (ii) of Prop. 1 accounts for the loss of one 
derivative in the inversion process. From Prop. 1, we can then obtain stability 
estimates of the form 

\\Tml\\H>'+^(n')+ X! II^P9')'llif''+i(o') + X! Il^op7llif'"(o') 

l<p<q<n-l l<p<n~l (18) 

The above stability estimate holds for fc = using the results of Proposition 1 and 
in fact for any fc > updating by standard methods (not detailed here [13]) the 
parametrices in (16) and (17) to inversions modulo operators in '^^^ (i.e., classi- 
cal VPDO of order — fc [13]) provided that the coefficients (70, {uj}) are sufficiently 
smooth. The presence of the constant C2 indicates that the reconstruction of 7 may 
be performed up to the existence of a finite dimensional kernel as an application of 
the Fredholm theory as in [16]. 

Equation (18) means that some components of 7 are reconstructed with a loss 
of one derivative while other components are reconstructed with no loss. The 
latter components are those that can be spanned by the components 2007 and 
{Tap^}i<a,(5<n-i- Somc algebra shows that the only such linear combination is 
Sr=o^ 7ii7, which, using the fact that X]"=o^ ^« ® ^ computed as 

^ T,,7 = (2^)-" / e"^-«Ao il„Ao 1 : 7(0 di 

- 7o~' : ((2^)-" e-«7(0 = tr (70"^), 

confirming the heuristics of Sec. 2.1. It can be shown that all other components 
of 7 (i.e. any part of 7^ in (8)) are, to some extent, spanned by the components 
Toa7, and as such cannot be reconstructed with better stability than the loss of one 
derivative in light of (18). Combining the above results with (14), we arrive at the 
main stability result of the paper: 

(7o'S)llHi-(f2') + ll7d||Hi-i(J2') ^ C'ljc'-H^ilffMO') + C'2il7llL2(n')- (19) 

Such an estimate holds for any fc > 1. 

The above estimate holds with C2 = when 7 dH{j) is an injective (linear) 
operator. Injectivity cannot be verified by microlocal arguments since all inversions 
are performed up to smoothing operators; see [3] in the isotropic setting. In the 
next section, we obtain an injectivity result, which allows us to set C2 = in the 
above expression. However, the above stability estimate (19) is essentially optimal. 
An optimal estimate, which follows from the above and the equations for (7, {vj}) 
is the following: 

l|A^lo7llffHf2') + hWH^-^in') < C\\dH\\Hk(n') +C2\h\\LHn') 

< C'(||A/|o7||h'=(o') + \h\\H^-Hn'))- 

The left-hand-side inequality is a direct consequence of (19) and the expression of 
dH. The right-hand side is a direct consequence of the expression of dH. The above 
estimate is clearly optimal. The operator M|o is of order 0. If it were elliptic, then 
7 would be reconstructed with no loss of derivative. However, M|q is not elliptic 
and the loss of ellipticity is precisely accounted for by the results in Lemma 2.3. As 
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we discussed above, it turns out that the only spatial coefhcient controlled by Af|o7 
is tr (7(7 "^7), and hence (19). 

2.3. Explicit inversion: Now, allowing 7 to be supported up to the boundary, 
we present a variation of the non-linear resolution technique used in [17, 18]. First 
considering n solutions generated by boundary conditions fulfilling Hypothesis 2.1, 
we establish an expression for 7 in terms of the remaining unknowns (ui, . . . , w„): 

7 = ^Q{[VU]H-^dHH-^[VUf - [\IV]H-^[VUf - [VU]H-^[VVf)-^Q, (20) 

where [VC/] and [VV] denote n x n matrices whose j-th columns arc Vuj and Vwj, 
respectively, and where H = {Hij}^^^^ and dH = In particular we 

find from (20) the relation 

tr ilo^l) ^ tr (H^^dH) - 2tr M, M {[VV][VUr^f . (21) 

Plugging (20) back into the second equation in (1) for 1 < i < one can deduce 
a gradient equation for the quantity tr (7(7^7) which in turn allows to reconstruct 
tr (7(7^7) in a Lipschitz-stable manner with respect to the LPD {dHij}'2j^i (i.e. 
without loss of derivative). 

Now turning to the full reconstruction of 7, we consider an additional solution 
Un+i generated by a boundary condition fulfilling Hyp. 2.2. The following propo- 
sition then establishes how to reconstruct {vi^ . . . , v„) from dH: 

Proposition 2. Assume that {gi, . . . , (7„+i) fulfill Hypotheses 2.1 and 2.2 over X 
and consider the linearized power densities dH = {dHij : l<i<j<n+l, i^ 
n + 1}. Then the solutions (wi, . . . , w„) satisfy a strongly coupled elliptic system of 
the form 

-V ■{-ioVv,) + W,j-Vv,^ f,{dH,V{dH)) (X), v,\ax ^ 0, I < i < n, (22) 

where the vector fields Wij are known and only depend on the behavior 0/70, Z and 
ui, . . . , u„, and where the functionals fi are linear in the data dHij. 

When the vector fields Wij are bounded, system (22) satisfies a Fredholm alter- 
native from which we deduce that if (22) with a trivial right-hand side admits no 
non-trivial solution, then (wi, . . . , w„) is uniquely reconstructed from (22). We can 
then reconstruct 7 from (20). 

Remark 3 (Case 70 constant). In the case where 70 is constant, choosing solutions 
as in Remark 1, one arrives at a system of the form (22) where Wij = if i 7^ j, so 
that the system is decoupled and clearly injective. 

The conclusive theorem for the explicit inversion is thus given by 

Theorem 2.4. Assume that {gi^ ... ,gn+i) fulfill Hypotheses 2.1 and 2.2 over X and 
consider the linearized power densities dH = {dHij : l<i<j<n-\-\^i^n + V\. 
Assume further that the system (22) with trivial right-hand sides has no non-trivial 
solution. Then 7 is uniquely determined by dH and we have the following stability 
estimate 

\\tr {1o^i)\\hHx) + I|7||l^(x) < C\\dH\\ H^(x)- (23) 
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2.4. Outline. We cover the microlocal inversion in Sec. 3. Linear algebraic and 
pseudo-differential preliminaries are given in Sec. 3.1. The leading-order symbols 
of order of the LPD functionals are computed in Sec. 3.2 and a proof of Lemma 
2.3 is given. The symbols of order —1 arc then computed in 3.3 and the proof 
Proposition 1 is given in Sec. 3.4. We then treat the explicit inversion in Sec. 4. 
Starting with some preliminaries in Sec. 4.1, we derive some crucial relations in 
Sections 4.2 and 4.3, before proving Proposition 2 and Theorem 2.4 in Sec. 4.4. 

3. Microlocal inversion. 
3.1. Preliminaries. 

Linear algebra. In the following, we consider the n x n matrices M„(R) with the 
inner product structure 

n 

A:B = tT (AB^) = ^^J^y-' (24) 

i,J = l 

for which M„(M) admits the orthogonal decomposition A„(R) © ^^(R). For two 
vectors U = (wi, . . . , Un)^ and V = {vi,. . . , Vn)^ in M" we denote hy U ®V the 
matrix with entries {uiVj^'^^^i, and we also define the symmetrized outer product 

U QV ■.^]^{U ®V + V ®U). (25) 
With • denoting the standard dotproduct on R", we have the following identities 
2UQV -.XQY ^{U ■X){V ■Y) + {U ■¥){¥ -X), U,V,X,Y (26) 
U ■ MU ^ M -.U ®U = M -.U QU, U e W\M e A/„(R). (27) 

Pseudo-differential calculus. Recall that we denote the set of symbols of order m on 
X by ^"(X), which is the space of functions p e C°°(X x R") such that for all 
multi-indices a and /3 and every compact set K C X there is a constant Ca,i3,K 
such that 

sup \dp,d^p{x,o\ < c^,^M^ + lei)"'"'"' 

We denote the operator p{x, D) as 

p{x,D)-t{x) = (27r)-" / e^--«p(x,m(C)d^ 

and the set of pseudo- differential operators (5*00) of order m on X by 
where 

^"^{X)^{p{x,D):peS"^{X)}. 

Suppose {?^^J■}o° is strictly decreasing and limmj = — oo, and suppose Pj e S"^''{X) 
for each j. We denote an asymptotic expansion of the symbol p e S''""(X) as 

p - ^ Pj G {X), for all k > 0. 

Given two VE'DO P and Q with respective symbols ap and erg and orders dp and 
dq, we will make repetitive use of the symbol expansion of the product operator 
QP = Q o P (see [11, Theorem (8.37)] for instance) 

aQp{x, - cjQUp + -V^aQ ■ V^ap + 0^"'^+"^-^), (28) 
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where 0(1(^1") denotes a symbol of order at most a. As we will need to compute 
products of three 'I'DO R. P and Q, we write the following formula for later use, 
obtained by iteration of (28) 

1 

<^RQP — <^Rf^Q<^P + -{^R^i^Q ■ ^x<7p + crgV^dfl; • V x<7p + O-pV^CTij • VxC^Q) 

(29) 

In the next derivations, some operators have matrix-valued principal symbols. How- 
ever we will only compose them with operators with scalar symbols, so that the 
above calculus remains valid. 

3.2. Symbol calculus for the LPD, properties of Mij\o and proof of Lemma 

2.3. Writing v^{x) = (27r)"" 4„ e'^-«w,(C) and 7(2;) = (27r)-" 4„ e'-^-«7(0 d^ 
(understood in the componentwise sense), we have 

Lov, := -V • (70 Vz;,) = (27r)-" / e'-< (^ • 70^ - t(V • 70) • C) MO d^, 

P,7 = V ■ (7Vu^) = (27r)-" J e^"-« (t^ Vu^ + V^u,) : 7(0 d^. 

Thus equation (4) reads LgVi = Pi-f, where the operators Lq := —V • (70V) and Pi 
have respective symbols 

(yLo=h + lu h := O 10^ and h := ■ -/o) ■ ^ e S\ (30) 

<JP. = P^,l + P^,o, K,i :=t?0VM. e (5^)"^" and k,o V^u, G (5")"^". 

(31) 

For Y a smooth vector field, we will also need in the sequel to express the operator 
y • V as 'I'DO, the symbol of which is denoted <jy v = fy.vli := • Y. 

We now write dHij as a '^DO of 7 with symbol Mij as in (9). dHij belongs to 
^>^(X) and we will compute in this paper the first two terms in the expansion of 
Mij (call them Mij \q and |_i), which in turn relies on constructing parametrices 
of Lq of increasing order and doing some computations on symbols of products 
of ^E'DO based on formula (28). If Q is a parametrix of Lq modulo 5'^'", i.e. 
K = QLq — Id € ^'~™, then straightforward computations based on the relation 
LoVi = Pi"f yield the following relation 

dHij{-i) = j:\/uiQ Vuj + (70 Vm^ • V) o Q o Pj^y + {joVuj • V) o Q o P,7 + Ky7, 

(32) 

where Kij := (70 Vu; • V) o KL^^Pj + {-faVuj ■ V) o KL^^P,. (33) 

For any i, Lq^ Pi denotes the operator 7 i-t- where Vi solves (4), and standard 
elliptic theory allows to claim that L^^Pi smoothes by one derivative so that the 
error operator Kij defined in (33) smoothes by m derivatives. In particular, upon 
computing a parametrix Q of Lq modulo the first three terms in (32) are 

enough to construct the principal part of the symbol Mij modulo 
Computation of My |o. In light of the last remark, we first compute a parametrix 
Q of Lq modulo 4'"^, that is, since Lq G Vl/^, we look for a principal symbol of the 
form ctq = (7_2 + 0(1^1^"^). Clearly, we easily obtain g_2 = ^2^^ = {£, ■ loO^^- 
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this case, the principal symbol of dHij at order zero is given by, according to (32) 
and (28), 

Mijlo = VUi VUj + (o-T-oV«.-v|l) 9-2 Pj,l + (o-7oV«j-v|l) 9-2 Pi,l 

= Vui Vuj - — ^— - ((7oVw, • Vuj) + (7oVuj • 0(^ © Vui)) ■ 
? • 70 ? 

Afij lo admits a somewhat more symmetric expression if pre- and post-multiplied by 
Aq^ the unique positive squareroot of 70, so that we may write. 



M,, 



, |o(x, ^ A-^ (v, Vj - iio ■ V,)io Vj - iio ■ V,)io 



A- 



(34) 



where we have defined ^0 ■= and a; := |x| for any x G M" — {0} as well as 
V, 



AoVui. This last expression motivates the proof of Lemma 2.3. 

Proof of Lemma 2.3. Proof of (i): Let 77 such that r/-^ = 0, and denote r/' 
so that ??' • Co = 0. Then using identity (26) and (34), we get 



2 \V, Vj - (lo ■ Vj - iio ■ V,)io V,\ : |o v' 
{V^ ■ ioWr V') + iV, ■ v'){V, ■ io) 

- do • v.)iio ■ ioWj ■ 11') - do • KO(c; • 7/)(y, • Co) 

- do • v,){i^ ■ ^^){v., ■ 77') - do • ^,)do • ??')do • V.,) 
0, 



where we have used ■^o • Co = 1 and Co • = 0, thus ('z^ holds. 
Proof of (ii): Recall that 



V, Vj - (Co • K)Co - (Co • VS )Co V, 
We write S',i(M) as the direct orthogonal sum of three spaces: 

5„(R) - (m Co Co) ({Co}"- {Co}^) © (Co {Co} 



^0 ^^^0 ^- 



(35) 



with respective dimensions 1, n(n— 1)/2 and n— 1. Decomposing A^^PA^^ uniquely 
into this sum, we write ^PtIq ^ ~ Pi + P2 + Pz- Direct calculations then show 
that 

M„-|o : ^0 '^^0 ' = ■■ i-Pi l<^<]< n. 

Since is a basis of E", {T^ Vj}i<i<j<„ is a basis of S'„(M) and thus (11) 

implies that 

-Pi + P2 = 0, i.e. Pi = P2 = 0. 

Therefore P = ^o^s^o with P3 = Co ??' for some 77' • Co = 0, so P = 70C r\ with 
77 proportional to A^t]' ^ i.e. such that 77 • C = 0, thus the proof is complete. □ 

In other words, all symbols of order zero A'/ij|o(a;, C) are orthogonal to the (x, C)- 
dependent n — 1-dimensional subspace of symmetric matrices 70C {C}"*"- One must 
thus compute the next term in the symbol exampansion of the operators dHij, i.e. 
Mijl^i. We will then show that enough symbols of the form Mijlo + Mij\_i will 
suffice to span the entire space «S'„(M) for every x G O' and C G S^, so that the 
corresponding family of operators is elliptic as a function of 7. 
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3.3. Computation of My |_i. As the previous section explained, the principal 
symbols Mij\o can never span S'„(M). Therefore, we compute the next term Mij\^i 
in their symbol expansion. We must first construct a parametrix Q of Lq modulo 
^I*^^, i.e. of the form 

aQ=q^2 + q-3 + Oi\^\-^), q,eS\ (36) 

Lemma 3.1. The symbols g_2 and g_3 defined in (36) have respective expressions 

g-2 - = (e ■ 7oO"\ (37) 
q~3 12^1^ Cp^gCi {[lo]pqdx, hohj - 2[johjdx^ bo]pq) ■ (38) 

Proof of Lemma 3.1. Using formula (28) with {Q,P) ^ {Q,Lo), and using the ex- 
pansions of aq and ctlq, wc get 

TQLo ~ + iq-2li + q-3h + ^V^g-a • ^^xh) + 0(1^1"^), 

In order to match the expansion 1 + + 0(1^1^^), the expansion above must satisfy, 
for large 

g_2^2 = 1 and q-2li + qsh + -'^^q-2 ■ "^xh = 0, 

L 

that is, q-2 = l2^ = ' loS,)~^ and 

q-3 = l2^ (^-q-2li - ^V5g_2 • "^xh^ = l2^ {-hh - ■ V xh) ■ 

Now, we easily have V4/2 = 270^ and Vxh = 9xj7o]pgCpCgei, where ei, . . . ,e„ is 
the natural basis of M". We thus deduce the expression of 53 

9-3 = '2^^'- {[lO\pq^p^qdxAlQ\ij^j - '^[lO\i]^]dxAlO\pq^p^q) : 

from which (38) holds, q-3 is clearly in from this expression, since l^^ is of 
order —6. The proof is complete □ 

We now give the expression of My|_i (or rather, that of ^0 ^^ij l-i ^o)- 

Proposition 3 (Expression of AoAfy|_iAo)- The symbol Mij\^i Aq admits the 
following expression for any {i,j) 

Ao M,,|_i(x,0 Ao = tllCoir' ((Co • Vj){M, - 2^0 H,^;) + Co H,F,) 

+ ^llCoir' ((Co • VA){mj - 2|o H.-Co) + io H,F,) 

+ l\M-' {v, ■ G(x,0(lo V.) + V. ■ Gix,0{io V,)) , 

(39) 

where we have defined Vi := A^Vui, Hi :~ Ao V^u^ A^, as well as the vector field 

G{x,0 ■■= llCof (t9-3Co + AoV,g_2) e {S"r. (40) 

Proof of Prop. (39). Assume Q is a parametrix of Lq modulo ^I*^^ and consider 
formula (32). Since the term 7 : Vui Vuj is of order zero, the computation of 
Mij I _ 1 consists in computing the second term in the symbol expansion oi RioQo Pj , 
and the same term with i, j permuted, where we denote Ri := 70 Vm^ • V with symbol 
ri,i = tjoVui ■ C. Plugging (t_r, = Tij, ctq = q^2 + q-3 and crp, = + Pi,o into 
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(29) and keeping only the terms that arc homogeneous of degree —1 in ^, we arrive 
at the expression 

crR,QPj\-i = rr,i{q-3Pj,i + q-2Pj,o) + -{pj,iV^ris ■ ^ xq-2 + V^{q-2ri^i) ■ "^xPjs)- 

(41) 

Note that the muhiphcations commute because the symbols of Q and Ri are scalar, 
while that of Pj is matrix- valued. Since My |_i = (^RiQPj |-i + o'RjQPi equation 
(39) will be proved when we show that 

Ao <yR.QP,\-i Ao - t|ieoir'((eo ■ v,)(m, - 2^0 q^L) 

. . . (42) 

+ Co O Hj-y, + V, ■ G{x, e)(eo Vj)j . 

Proof of (42). Starting from (41), plugging the expression r^^i = t{Vi-^o), using 
the identity 

Vj((?_2?'i,l) • "^xPj,! = O (V^WjV5((7_2?'i,l)), 

and pre- and post-multiplying by Aq yields the relation 

^0 crR.QP,\-i Ao = l{V^ ■ £.o){q-3L^o & Vj + g^aHj) 

(43) 

+ {V, ■ AoV:,g_2)tCo Vj + Co HjA„ 1 V4(g_2r,4). 

Gathering the first and third terms recombines into t||Co||~^^i • G(Co Vj) (the last 
term of (42)). On to the second and fourth terms, we first compute 

ivar,i'?-2) = t(^. • Co)(-iieoir')2eo + iiCoir'tK: = 4ur\v, ~ 2{v, ■ eo)io). 

Using this calculation, the second and fourth terms in (43) recombine into 

^lieoir' ((^0 • K:)(Hj - 2|o H,lo) + HjK:^ 
thus the argument is complete. □ 

3.4. Proof of Proposition 1. Preliminaries: By virtue of Hypothesis 2.1, Vu„+i 
may be decomposed into the basis Vui, . . . , Vu„ by means of scalars /ii, . . . , /i„, /i 
such that 

n 

y ^Vu, + Vu„+i = 0. (44) 

As seen in [17, 18], the coefficients fii, . . . ,/in+i are directly computible from the 
power densities {Vui-7oVMj}i<i<j<„+i and on the other hand, we have the relation 

i 

fii det(Vui, . . . , Vu„+i, . . . , Vw„) ^ . ^ 
fj. det(Vui, . . . , Vum) 

thus ~ Zi as defined in (7) for I < i < n. In the next proofs, we will use the 
following 

Lemma 3.2. Under hypotheses 2.1 and 2.2, the following matrix-valued function 

M := i^iMt + (45) 
is symmetric and uniformly invertible. 
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Proof. Symmetry of M is obvious by definition. Taking gradient of (44) , we arrive 
at 

n 



Pre- and post-multiplying by Aq, we deduce that 

M = fi,m, + AiH„+i = -fiAnZ, ® = -fj.AnZV'^, 

where V := . . . \ Vn\- The proof is complete since Hyp. 2.1 ensures that fi never 
vanishes and V is uniformly invertible, and Hyp. 2.2 ensures that Z is uniformly 
invertible. □ 

The Tpq operators. Proof of Prop. 1: As advertised in Sec. 2.2, because 
of the algebraic form of the symbols of the linearized power density operators, it 
is convenient for inversion purposes to define the microlocal change of basis Tj = 
{Tpqj}i<p<q<n as lu (13), i.c. 

Tpq^ix) := (27r)-" / e'-^'^A^'ip iqA^' : m d^- 

To convince ourselves that this collection forms a microlocally invertible operator 
of 7, let us introduce the zero-th order ^'DOs Pijpq with scalar principal symbol 
o'-Pijp, ■= (^i ■ AoS.p)iej ■ AoCq) for 1 < i,j,p,q < n. Then for any 1 < i < j < n, the 
composition of operators Pijpq ° Tpq has principal symbol (repeated indices 

are summed over) 

(e, ■ Aoip){ej ■ Aoiq)A^^^p ^ = (e, ■ ^o4)^o ^4 © i^j ' ^o4)^o ^4 

where we have used the following property, true for any smooth vector field V: 

V^iV- AoQA-'ip. 

Thus for any I < i < j < n, the composition X)p g=i Pijpq ° Ppq recovers 7^ = 7 : 
e^ ej up to a regularization term. This in particular justifies the estimates (14) 
and the subsequent inversion procedure. We are now ready to prove Proposition 1. 

Proof of Proposition 1. From the fact that (Vi, . . . , X^) is a basis at every point and 
given their dotproducts Hij = Vi ■ Vj , we have the following formula, true for every 
vector field W: 

W = H^^W ■ Vp)Vq. (46) 

Proof of (i): Reconstruction of the components T007 and {Tap^}i<a<i3<n-i- 
We work with Afylo Aq^V^Q V, ~ (^o • V,% V,- - (^^0 • V,)io V,. 

Using (46) with W = straightforward computations yield 

• ^9)^"" (6 • VpMj lo = L ip - iio ■ L)io 6 - (^"o • ip)io L 

i,3,P,q 

' -Co 0^0 if a = /3 = 0, 
if = Q' ^ /3, 
L^ip if a^0,/^^0. 
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which means that upon defining Qapij G with scalar principal symbols 

'^Qoo., ■.= ~Y.H''{io ■ V,)HP^io ■ Vp), 
p,q 

^Q.,., E ^" (^" • -Vpl 1 < « < /3 < n - 1, 

p,q 

relation (16) is satisfied in the sense of operators since the previous calculation 
amounts to computing the principal symbol of the composition of operators in (16). 
Proof of (ii): Reconstruction of the components {7oq7}]^<q,<„_j^. It remains 

to construct appropriate operators that will map dH{'-f) to the components Toaj for 
1 < a < n — 1, which is where the additional measurements dHi n+i come into play. 
Let (^1, . . . ,fJ.n,f^) as in (44) and construct the Vl/DO {^4(7)}"=! as in (15). It is 
easy to see that, since the /x, are only functions of x, the terms of fixed homogeneity 
in the symbol expansion of L; satisfy 

(TL,\k = fJ-jMijlk + ^Afj,„+i|fe, fc = 0, -1, -2, . . . . 

Then from equation (34) and relation (44), we deduce that ctl; |o = 0, so that 
Li G ^P^^. Moreover, using equation (39) together with relation (44), we deduce 
that 

ffi. |_i ^ Ao CJL, 1-1 = illCoir' ((Co ■ ^^)(M - 2^"o M^o) + ^"0 © MV^,) 

is now the principal symbol of Li. Using relation (46) with W = M^^^q,, the 
symmetry of M and multiplying by M, wc have the relation 

Using this relation, we deduce the following calculation, for 1 < a < n — 1 

HP\io. ■ M-^Vp) dL, 1-1 = tllCoir' (do • M-iea)(M - 2|o M|o) + |o L) ■ 

(47) 

While the second term gives us the missing components ToqT, we claim that the 
first one is spanned by ^0 $0 and {^q ip}i<a<(}<n-i- Indeed we have 

(M - 2^0 Me'o) : (Co e'a) - 0, 1 < a < n - 1, 
(M - 2|o M|o) : do Co) = -|o • M|o, 

(M - 2^'o M^^o) : {L 6) = L ■ M^^, l<a</3<n-l, 
so we deduce that 

M - 2^0 M^'o = - iio ■ Mio ) Co + E ' ^^"^ ) (4^) 

l<a,P<n-l 

In light of these algebraic calculations, we now build the parametrices. Let 

i| e ^-l, e ^-O, i? e (^rO)«xn^ g ^0 j^j^ji g ^0 ^^j, ^J30s with 

respective principal symbols 

a 1 = -l\M, tTB„ = HP'iia. ■ M- Vp), aR^M- 2^0 M|o, 
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Then the relation (47) imphes (17) at the principal symbol level. The operator R 
can indeed be expressed as the following zero-th order linear combination of the 
components Too and {Tai3}i<a,p<n-i- 

R = —RooTqo + ^ RapTap 
l<a,P<n-l 

'RooQooij + ^ RajsQapij o dH^j mod^*"^, 

l<Q.,8<n-l j 

SO that the left-hand side of (17) is expressed as a post-processing of measurement 
operators dHij only. The proof is complete. □ 

4. Explicit inversion. 

4.1. Preliminaries and notation. For a matrix A with columns Ai, . . . , An and 

(ei, . . . , e„) the canonical basis, one has the following representation 

n n 

A = ^^Aj'gj Bj and = ^ e^ Aj. 

i=i j=i 

More generally, for two matrices A = [Ai \ . . . and B = [Bi \ . . . |i?„], we have 
the relation 

n 

^Aj(gBj = AB^. 
j=i 

Finally, for A a matrix and V ^ [Vi \ . . . \ Vn\, the sum AijVj is nothing but the z-th 
column of the matrix VA^. 

4.2. Derivation of (20) from Hypothesis 2.1: Let us start from n solutions 
(mi, . . . , Un) fulfilling Hypothesis 2.1, and let (wi, . . . , v„) the corresponding solutions 
of (4). We also denote [VU] := [Vui| . . . |Vu„] and [W] similarly. We first mention 
that for any vector field V, we have the following formulas 

V = HPi{V ■ 7oVMp)Vug = HP'^iV ■ Vup)7oVug, (49) 

which also amounts to the following matrix relations 

HP'^iVup (g) Vug)-fQ = RP'^-foiVup (E) Vug) = I„. (50) 

From the relation 

dH^j = (7Vui -I- 7o Vwi) • Vuj + 70 Vu^ • Vu^, I <i,j <n, 

we deduce, using (49), 

7VU, + 7oVwi = HP'^ (dHip - 7oVwp ■ Vu,) 70 Vw,, 1 < i < n. (51) 

The previous equation allows us express 7 in terms of the remaining unknowns 
(wi, . . . , Vn)- Indeed, taking the tensor product of (51) with H^^joVuj and summing 
over i yields 

7 + 7oVwi (g) yuj-foW^ = HP1 {dH,p - 7oVi;p • Vu;) (70 Vu, ® VujjoH'^) 

= dH,p-/o{HP''yuq ® i/*^Vuj)7o - 70 Vu, <g) Vvp^/oHP", 
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where we have used the identity (49) in the last right-hand side. We may rewrite 
this as 

7 = 70 {dH,p{HP'^yuq ® H'^Vuj) - 2W^Vvi Vuj) 70. (52) 

One may notice that the above expression is indeed a symmetric matrix. In matrix 
notation, using the prehminaries, we arrive at the expression (20). 

4.3. Algebraic equations obtained by considering additional solutions: 

Let us now add another solution Un+i with corresponding solution Vn+i at or- 
der 0{e). By virtue of Hypothesis 2.1, as in section 3.4, Vm„+i may be expressed 
in the basis (Vui, . . . , Vu„) as 

n 

V— Vu, -I- Vu„+i =0, (53) 

ti ^' 

where the coefficients fii can be expressed as ratios of determinants, or equivalently, 
computable from the power densities at order e", see [17, Appendix A. 3]. For 
I < i < n, we define Zi V(/i~^/ii), and notice that we have the following two 
algebraic relations 

n n 

Z, ■ 70 Vu., = and ^ Z^ A dui = 0. (54) 

i=l 1=1 

The first one is obtained obtained after applying the operator V • (70') to (53) and 
the second one is obtained after applying an exterior derivative to (53) . 
Moving on to the study of the corresponding Vn+i solution, we write 

dHn+ij + —dH,,j = f Vi;„+i + —Vu, J ■ 7oVMj, 1 < j < n, 
V / 

where we have cancelled sums of the form (53). Using the identity (49), we deduce 
that 

V-y^+i + (/i- Vi)Vw» = (di/„+i,p + {fi-^fi,)dHip) Vug. (55) 
Taking exterior derivative of the previous relation yields 

Z^^ Adv^ = d {HP''{dHn+i,p + {fi-^n^)dH^p)) A dug. (56) 
We now apply V • (70-) to (55), the left-hand side becomes 
V • (7o(Vu„+i + (/.t"Vi)Vwj)) . . . 

= V • (7oVw„+i) -I- Zi ■ 7oVwi + {jjL^^iii)V ■ (70 Vw,;) 
= -V • (7Vw„+i) -I- Zi ■ -faVvi - ■ {'^Vui) 

= Zi ■ 7oVui - V • (7(Vu„+i -I- {^i^^ ^jLi)Vui)) + Zi ■ jVu^ 
= • (70 Vu; + 7VMj), 
thus we arrive at the equation 

Z^ ■ (70 Vw, + jWu,) ^ V {HPi{dHn+up + {^1-^ fii)dH,p)) ■ 70 Vu, Yg ■ 7oVu„ 
where the vector fields 

Yg V {HP''{dH^+i,p + Ur^fi,)dH,p)) , 1 < g < n, (57) 
are known from the data dH. Combining the latter equation with (51), we obtain 
{Zi ■ ^QVug)HP'^ {dHip - 7oVwp ■ Vui) = Yg ■ 7oVMg, 
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which we recast as 

The left-hand side can be considerably simplified by noticing that the second equa- 
tion of (54) implies [VJ7]Z^ = Z[VJ7]'^. With this fact in mind, the left-hand side 
looks like Xp ■ Vup, where wc compute 

Xp = HP'^-foVu, ® Zao^ug = 7o[VC/]Z^7o[VC/]i/-iep 

= joZ[VUf[\/U]-^ep = 7oZp. 

Finally, wc obtain the more compact equation 

n 

Y,loZp-Vvp^ f, where f -.^ {RP^dH^p Z, ~ Y^) ■ ^o'^Ug, (58) 
p=i 

with Yq given in (57). 

Remark 4 (On algebraic inversion). In equations (56) and (58), the only unknown 
is the matrix [WV] := [Vwi, . . . , Vw„]. Equations (56) and (58) give us the projection 
of that matrix onto the space Z^„(R) and onto the line WjoZ respectively. As in 
the non-linear case [17, 18], we expect that a rich enough set of such equations 
provided by a certain number of additional solutions (u„_|-i, . . . ,ii„_|-/) leads to a 
pointwise, algebraic reconstruction of [Vy], however wc do not follow that route 
here. 

4.4. Proof of Proposition 2 and Theorem 2.4. We now show that provided 
that we use one additional solution (on top of the basis (ui,...,u„)) such 

that the matrix Z is of full rank, then we can reconstruct (wi, . . . , w„) via a strongly 
coupled elliptic system of the form (22), after which we can reconstruct 7 from 
{Vvi, . . . , Vvn) by formula (20). We now show how to derive this elliptic system. 

Proof of Proposition 2. According to Hypothesis 2.2, the matrix Z = [Zi \ . . . |Z„] 
has full rank and we recall the important equations 

n n 

7oZp ■ Vwp = / and Z\ A dvi = cj, where (59) 

p=l i=l 

LJ = Y^A dug, Yq := V{HP^idHn+i,p + ^i,)dH,p)), (60) 

and where / is given in (58). Assuming that Z has full rank, the family (Zi, . . . , Z„) 
is a frame with dotproducts defined as S^- = Zi ■ Zj^ and in this case we define its 
dual frame Z* := for 1 < i < n, such that Z* ■ Zj = Sij, i.e. with Z* the 

matrix with columns Z* , we have the relation Z* — Z^^ . The second equation of 
(59) may be rewritten as 

ZI ■ Vvp - Z; ■ Vvq = L^(Z;, Z*), l < p, g < n. (61) 

Applying the differential operator Z* ■ V to the first equation of (59), we obtain 

n 

Y^{Zl ■ V)(7o^p • V)«p = {Z: ■ V)/. (62) 
p=i 
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Using (61), we may rewrite the left-hand side of (62) as 

(Z* ■ V)(7oZp • V)vp = [Z*, 7oZp] • Vvp + i-foZp ■ V)(Z* • V)vp 

= [z*, jnZp] ■ + {-foZp ■ v){z; ■ V)v, . . . 

+ (7o^p-V)(c.(Z;,Z*)), 

where we have introduced the Lie bracket of two vector fields, which may be written 
in the Euclidean connection 

[X, Y] := {X ■ V)Y - {¥ ■ V)X. (63) 

Plugging the last calculation into (62) (repeated indices are summed over) 

(7o^P • V)(Z; • V)^;, + [Z*, -foZp] ■ Vvp = {Z* ■ V)/ - (70^^ ■ v)(c.(z;, Z*)). 

(64) 

We now look more closely at the principal part of this equation. The first term may 
be written as 

7oZp ® z; : \/\, + ii^oZp ■ V)z;) ■ Vf, = 70 : v^v, + {{joZp ■ v)z;) ■ \/v„ 

where we have used that Zp (g) Z* ~ I„ . We thus obtain a strongly coupled elliptic 
system of the form (22), where 

W,, (V • 70 - iiloZp ■ V)Z;)) 5,, - [Z*,7oZ,], 1 < z, j < n, (65) 

h~~Z*-Vf+{-ioZp-V){uj{Z;,Z*)), l<i<n. (66) 

This concludes the proof. □ 

In order to assess the properties of system (22), we recast it as an integral equa- 
tion as follows: Let us call Lq := —V ■ (70V), and define L^^ : H^^{X) 3 f u E 
Hq{X), where u is the unique solution to the equation 

-V-(7oV?.)-/ (X), u\ax=0- (67) 

By the Lax-Milgram theorem (see e.g. [10]), one can establish that such solutions 
satisfy an estimate of the form ^ C'll./'l|ff-i(x): where C only depends on 

X and the constant of ellipticity of 70, thus Lq^ : H^^{X) — > H^{X) is continuous, 
and by Rellich imbedding (i.e. the fact that the injection H^^ is compact), 

I/Q ^ : L'^{X) H^{X) is compact. 

Applying the operator Lq^ to (22), we arrive at the integral system 

n 

^^ + Yl ^o\W,j ■ \/vj) = /i, := V, {X), l<i<n, (68) 

where it is easy to establish that for I < i,j < n, the operator 

: H^iX) 3v^ P,^v := L^\W,, ■ Vv) e H^{X) (69) 

is compact whenever the vector fields Wij are bounded. In vector notation, if we 
define the vector space H = (i/g (X))", v = (wi, . . . , u„), h = {hi, . . . , /i„) and for 
V e ^, 

Pv := (Py w„ F2jVj, • ■ • , PnjVj) e n, (70) 

we have that P : H ^ 7^ is a compact linear operator, and the system (22) is 
reduced to the following Fredholm (integral) equation 

(I + P)v = h. (71) 
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Note here that the operator P defined in (70) depends only on 70 and the solu- 
tions Mi, so that the injectivity properties depend on the 70 around which we pose 
the problem, in particular, whether one can fulfill hypotheses 2.1 and 2.2. 
Injectivity and stability. Equation (71) satisfies a Fredholm alternative. In particu- 
lar, if —1 is not an eigenvalue of P, (71) admits a unique solution v G H (injectivity), 
(I + P)^^ : H — > H is well-defined and continuous and v satisfies the estimate 

||v||«< ||(I + P)-i||,:(«)|lh||„, (72) 

from which we deduce stability below. In the statement of Theorem 2.4, the fact that 
"system (22) with trivial right-hand sides admits no non-trivial solution" precisely 
means that —1 is not an eigenvalue of the operator P. 

Remark 5 (Injectivity when 70 is constant). When 70 is constant, constructing 
{ui, . . . , Un+i) as in Remark 1 yields Z = Q a constant matrix. In particular, the 
commutators [Z*,jQZj] vanish in the expression (65) of Wij. Thus system (22) is 
decoupled and clearly injective. By continuity, we also obtain that (22) is injective 
for 7o (not necessarily scalar) sufficiently close to a constant. 

We now prove Theorem 2.4. 

Proof of Theorem 2.4- Starting from the integral version (71) of the elliptic system 
(22) in the case where —1 ^ sp(P), then the Fredholm alternative implies (72). In 
order to translate inequality (72) into a stability statement, we must bound h in 
terms of the measurements {dHij}. We have for 1 < i < n, 

\\hi\\H^{X) < \\Ln^\\c(,H-\Hl,)\\fi\\H-^X), 

and since /,;, expressed in (66) involves the dHij and their derivatives up to second 
order, if we assume all other multiplicative coefficients to be uniformly bounded, 
we obtain an estimate of the form 

\\h,\\Hi^x) < C\\dH\\Hi(x), where \\dH\\H^x) ■= \\dH,j\\HHx), 

thus wc obtain in the end, an estimate of the form 

MHiix)<C\\dH\\H^x)- (73) 

Once V is reconstructed, we can reconstruct 7 uniquely from dH and [VV] using 
formula (20), with the stability estimate 

ML^X)<C\\dH\\H^^x)■ (74) 

Regaining one derivative back on tr (7(7^7): In order to see that tr (7(7^7) 
satisfies a gradient equation that improves the stability of its reconstruction, the 
quickest way is to linearize [18, Equation (7)] derived in the non-linear case, which 
reads as follows: 

Vlogdet7- = Vlogdeti/'^ + 2 ((V(i?'')^') ■7''Vuf) Vu^, 

where is the n x n matrix of power densities Hfj = \7uf ■ 7^Vu^ and {H'^y^ is 
the (j, l)-ih entry of (H^)^^. Plugging the expansions 7"^ = 70 + £7, uf = Ui + evi, 
Hfj = Hij + edHij , and using the fact that 
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the linearized equation at 0(e) reads 

ivtr (7o- V) = ^Vtr (H-^dH) + [VW' ■ 7oVu/)Vw, + [VW' ■ 7oVfz)VMj 

From this equation, and using the stabihty estimates (73) and (74), it is straighfor- 
ward to cstabhsh the estimate 

l|tr ho'l)\\HHX)<C\\dH\\ 

and thus the proof is complete. □ 
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